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Abstract. - Gyration radius of individual's trajectory plays a key role in quantifying human 
mobility patterns. Of particular interests, empirical analyses suggest that the growth of gyration 
radius is slow versus time except the very early stage and may eventually arrive to a steady 
value. However, up to now, the underlying mechanism leading to such a possibly steady value 
has not been well understood. In this Letter, we propose a simplified human mobility model 
to simulate individual's daily travel with three sequential activities: commuting to workplace, 
going to do leisure activities and returning home. With the assumption that individual has 
constant travel speed and inferior limit of time at home and work, we prove that the daily moving 
area of an individual is an ellipse, and finally get an exact solution of the gyration radius. The 
analytical solution well captures the empirical observation reported in [M. C. Gonzalez el al., 
Nature, 453 (2008) 779]. We also find that, in spite of the heterogeneous displacement distribution 
in the population level, individuals in our model have characteristic displacements, indicating a 
completely different mechanism to the one proposed by Song el al. [Nat. Phys. 6 (2010) 818]. 
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Understanding human mobility patterns is important to 
many different disciplines from epidemiology to communi- 
cation engineering and transportation engineering [IH1]- 
With the development of modern electronic techniques, 
some novel devices including online bill tracker [5], mo- 
bile phone [HIE] and Global Positioning System [5J|2] have 
been used to track human mobility trajectories. By an- 
alyzing these data, scientists have revealed some remark- 
able features of human mobility patterns, such as heavy- 
tail distributions of travel distances and waiting times, 
spatio-temporal regularity and bounded nature of indi- 
vidual trajectory [5till| . To explain the origin of observed 
scaling properties in human mobility patterns, formula- 
tion of simplified models and their quantitative analysis 
have attracted increasing interests from many branches of 
sciences [T2HT7] . 

A key quantity related to human mobility trajectories is 
gyration radius |18j , a measure of how far from the center 
of mass the mass is. Gonzalez, Hidalgo and Barabasi [5] 
used gyration radius to quantify individual trajectory 



tracked from numerous anonymized mobile phone users. 
The gyration radius grows fast in the very early stage and 
then become very slow versus time. It approximately ap- 
proaches to a steady value in hundreds of hours. They sug- 
gested a logarithmical function to fit the empirical data. 
Moreover, they found that gyration radius has a strong 
impact on travel distance distributions over all users. Af- 
ter the removal of the dependence of gyration radius, the 
travel distance distributions of groups collapsed into a sin- 
gle curve. Therefore, they suggested using gyration radius 
as a characteristic travel distance for each individual. In 
a word, these empirical evidences indicated that the gy- 
ration radius of individual's trajectory plays a key role in 
characterizing human mobility patterns. 

Despite its importance for quantitatively understanding 
human mobility patterns, the factors responsible for the 
gyration radius of trajectory is still unclear. Several mod- 
els Q3HH] have been developed to reproduce the spatio- 
temporal patterns of human mobility, but none of them 
explained the mechanism stabilizing gyration radius. Re- 
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cently, Song et al. [12] proposed a novel model for individ- 
ual human mobility. Their model is based on two compet- 
itive mechanisms: exploring unknown locations and re- 
turning preferentially to familiar places. The model ex- 
plained many scaling laws observed in human mobility 
patterns, as well as the ultraslow increasing phenomenon 
of gyration radius. However, the factors affecting the gy- 
ration radius of individual's trajectory have not been fully 
documented. 

In this Letter, we proposed a simple model for human 
mobility as well as a method to solve it. This model is 
motivated by empirical evidence drawn from a travel di- 
ary data set named Mobidrive [!§]• The data set records 
the day-to-day travel behaviors of 360 persons in two Ger- 
many cities over a six-week survey period. A series of 
researches based on this data set have revealed some com- 
mon features of individuals' daily travel behaviors [2Url2"3"] . 
The most notable one of these features is the similarity of 
individual's weekday travel activities [2TJ1I2T] , which means 
that most individuals take their working-day activities in 
nearly the same order: commuting to workplace (or school 
for students) in the morning, spending their daytime at 
workplace, doing some leisure activities (shopping, dining, 
entertainment, etc.) and returning home in the evening. 
Although the times and visiting locations of leisure ac- 
tivities can be variable in this "home-work-leisure-home" 
trips chain [35], the activities that periodic commuting to 
workplace and returning home are nearly invariant. Sim- 
ilar regular patterns were also observed in some mobile 
phone users' travel behaviors fM\ . 

Here we introduce a straightforward model to reproduce 
the regular travel behavior of human individual in working 
days. For simplicity, we assume that: (i) each individual 
has fixed locations of home and workplace; (ii) each in- 
dividual has an inferior limit of time at home and work 
everyday; (iii) each individual has a constant travel speed 
relied on her/his common means of transportation; (iv) 
each individual takes one trip per day(A trip specifically 
stands for the travel from workplace to a selected loca- 
tion for leisure activities. The other two movements, from 
home to workplace and from the selected location to home, 
are not accounted). The last assumption is based on the 
result of a previous empirical research on Mobidrive |22j . 
which indicated that in most time an individual performs 
no more than one trip per day. Based on the assumptions, 
we present a detailed description of our model: in every 
working day, an individual first commutes to workplace 
from home, then goes to a random location for leisure 
activities, and finally returns home. The individual can 
select any locations to do after- work activities, but she/he 
must guarantee to return home in time. In other words, 
because of the restrictions of travel speed and the maxi- 
mum possible leisure time everyday, the individual is pro- 
hibited to select the points far away from home and work- 
place in our model. Then, after T simulation steps, we will 
obtain a trajectory composed with T points at "home" , T 
points at "workplace" and T points scattered in a defined 
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Fig. 1: (Color online) Elliptical area of an individual's daily 
travel. Two foci of the ellipse are home and workplace, and 
the distance between them are focal length (2c). The major 
axis length (2a) represents the maximum daily travel distance. 



area. 

The defined area of individual's trajectory in our model 
can be shaped by a geometric approach: with the model 
rules described above, the maximum possible travel dis- 
tance in one's everyday after-work activities (including 
home-return) is 



2a = v(24 - t rn ) - 2c, 



(1) 



where v is the travel speed, t m is individual's inferior limit 
of time at home and work everyday, c is half of the dis- 
tance between home and work (c < a). The time resolu- 
tion is hour. Obviously, we can draw the maximum area 
of individual daily travel as an ellipse (see Fig. [IJ . The 
ellipse's two foci are individual's home and workplace, the 
focal length is the distance between home and workplace, 
and the major axis length, 2a, is the maximum possible 
daily travel distance subtracting commute distance. The 
elliptical area of individual human mobility has been em- 
pirically found in Mobidrive [53] (see also Fig. 3a in Ref. [5] 
where the mobile phone users' trajectories seem to be also 
limited in elliptical areas). 

With the obtained elliptical area, we next solve the gy- 
ration radius of trajectory for our model. The gyration 
radius of the trajectory can be expressed as 
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where T is the simulation step, r cm is the center of mass 
of the trajectory, and in each time step, three locations 
are considered, namely r i = r 4 = r 7 = • • • = home, r*2 = 
f*5 = rg = ■ • • = workplace, and r3,re,rQ, • ■ ■ , r^T are the 
locations randomly selected in the ellipse. When T is large 
enough, according to the geometrical symmetry, we can 
deduce that the center of mass of the locations is at the 
center of the ellipse. Notice that the 2T points at foci 
(home and work) of the ellipse have same distance, c, to 
the center of the ellipse, we can rewrite Eq. to 
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Fig. 2: (Color online) The time evolution of gyration radius in 
the model (c = 30). The curves are numerical simulations for 
different a, and the straight lines represent analytical results. 
Data points are averaged by 100 independent runs, each lasts 
200 steps after initialization. 



where I is the mean square distance from T random points 
to the center of the ellipse. This distance can be calculated 

as 



_ _ AO J^(a 2 cos 2 6 + b 2 sm 2 9)abp 3 dp _ a 2 + b 2 
■nab 4 



(4) 



where b is the half minor axis length of the ellipse. Com- 
bining Eqs. ([2]|4]) and b 2 = a 2 — c 2 , we finally obtain 



7c 2 + 2a 2 
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(5) 



The result shows that the gyration radius of individ- 
ual's trajectory relates to two factors: one is the distance 
between individual's home and workplace, another is in- 
dividual's maximum daily travel distance, which depends 
on individual's travel speed and daily time assignment. 
The result also implies that the individual who has more 
leisure time and more advanced means of transportation 
will has a larger area of daily mobility. 

Notice that our analytical result shows that the radius 
of gyration is a constant, but the previous empirical obser- 
vation [5] shown that it grows with time in the very early 
stage. To explain this difference, we simulate the time 
evolution of gyration radius for our model (see Fig. [2]) 
. The simulation results demonstrate that the gyration 
radii grow very fast in the early stage and will converge 
asymptotically to stable values in accordance with the cor- 
responding solutions. Although one may choose logarith- 
mic function to fit the simulation curves before they con- 
verged, the curves are indeed in the relaxed stage. The 
empirical observations reported in Ref. [6] look like the 
simulations shown in Fig. [2] Although these empirical ob- 
servations can be well captured by logarithmic functions, 
they may converge to a steady value after a very long time, 
as suggested by the present model. In addition, it is more 



Fig. 3: (Color online) Travel distance distributions of three 
individuals with the same r g = 15.12 in the model. Data points 
are averaged by 100 dependent runs, each run lasts 1000 steps 
after initialization. 



convenient to use a single value rather than the whole of 
the time evolution curve of gyration radius to characterize 
human mobility patterns. 

Additionally, from Eq. (|5|) we know that individuals 
sharing same r g can have different c. The parameter c 
have crucial influence on travel distance distribution of 
individual because there are one third of trips taken be- 
tween home and workplace in our model. To demonstrate 
this, we simulate the travel distance distribution of indi- 
viduals with same r g but different c, and plot the results 
in Fig. [3] As the figure shows, all of the distributions 
arc different: they peak respectively at 2c, and, none of 
them is right-skewed. This seems incompatible with the 
empirical observation that human have heavy-tail distri- 
butions of travel distances [5,6 . Nevertheless, to the best 
of our knowledge, there were not sufficient evidences show- 
ing that the travel distance distribution of individual was 
heavy-tailed. In spite of the empirical evidence shown 
that the population with same r g has heavy-tailed dis- 
tribution of travel distances, it is insufficient to deduce 
the existence of a similar scaling law for individual. In 
the consideration of only daily travel instead of long-term 
traveling by trains, airplane, etc, • ■ •, pealed law at indi- 
vidual level was also be observed [25] . Actually, the emer- 
gence of heavy-tailed nature at the population level can 
also be resulted from the heterogeneity of the distance be- 
tween home and workplace (such heterogeneity has been 
reported in Ref. [HH?]). 

Indeed, the empirical analyses at the individual level 
[2"5] suggested that most individuals' movements are lim- 
ited in an ellipse-like area with a peaked displacement dis- 
tribution. Figure 4 reports an example where the dis- 
tributions of c follows -P(c) <~ c _1 , where the displace- 
ment distribution can be approximately fitted with double 
power-law form. The result indicates that the heteroge- 
neous displacement distribution emerged at the popula- 
tion level can be resulted from the heterogeneity of home- 
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Fig. 4: (Color online) Travel distance distribution of population 
with the same r g — 200 but different c following P(c) ~ c _1 in 
the model. Data points are averaged by 100 independent runs, 
each lasts 1000 steps after initialization. 



workplace distances or some other factors instead of the 
heterogeneous displacement distribution at the individual 
level. 

In conclusion, we have proposed a simplified human mo- 
bility model and given an exact solution of gyration radius 
for it. The analytical solution agrees well with the empir- 
ical observation of gyration radius of mobile phone users, 
suggesting that our model may correctly capture the main 
mechanism leading to the bounded nature of human mo- 
bility. Different from the random diffusion of particles, 
human daily travel has a high degree of periodicity and 
regularity, which is the root of the high predictability of 
individual movements and the rhythms of urban traf- 
fic Additionally, individuals in urban travel are pur- 
poseful, or socially contextualized, and they have differ- 
ent destination choice strategies corresponding to various 
trips purposes. Our model, although simplified, incorpo- 
rates these intrinsic periodicity and purposivencss of hu- 
man travel behavior, and is expected to be extended to 
a more realistic traffic prediction approach in future re- 
search. 
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